Effects of density imbalance on the BCS-BEC crossover in semiconductor 

electron-hole bilayers 
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We study the occurrence of excitonic superfluidity in electron-hole bilayers at zero temperature. 
We not only identify the crossover in the phase diagram from the BCS limit of overlapping pairs 
to the BEC limit of non-overlapping tightly-bound pairs but also, by varying the electron and hole 
densities independently, we can analyze a number of phases that occur mainly in the crossover region. 
With different electron and hole effective masses, the phase diagram is asymmetric with respect to 
excess electron or hole densities. We propose as the criterion for the onset of superfluidity, the jump 
of the electron and hole chemical potentials when their densities cross. 
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PACS numbers: 03.75.Hh,03.75.Ss,73.63.Hs 

There is a growing interest in the physics of the BCS- 
BEC crossover, owing to the recent experimental ad- 
vances with trapped Fermi atoms. With the use of Fano- 
Feshbach resonances, this crossover has been observed 
with fermionic 6 Li and 40 K atoms |l| which become com- 
posite bosons in the BEC limit. Although trapped atoms 
represent an ideal testing ground for a fundamental un- 
derstanding of the BCS-BEC crossover, technological ap- 
plications exploiting the occurrence of condensates will 
most probably rely on semiconductor systems. In these 
systems, excitons made up of electrons and holes play the 
role of composite bosons. 

Excitonic systems were, in fact, the first to be consid- 
ered for the BCS-BEC crossover. In bulk materials, orig- 
inal work on exciton condensation was done by Keldysh 
and co-workers Q. Extension of this to the BCS-BEC 
crossover was proposed by Nozieres and Comte Q . How- 
ever, in bulk materials fast electron-hole recombination 
hinders the detection of BEC. Proposals have accordingly 
been made to condense excitons with spatially separated 
electrons and holes 0, . Particularly promising are bi- 
layer quantum-well systems separated by a distance d, 
with conduction-band electrons in one well and valence- 
band holes in the adjacent well |E IS El • Recently, de- 
velopments for the detection of excitonic BEC in coupled 
quantum wells has been reported [TnL fill Hp. Techno- 
logical advances with electron-hole bilayers [13j make it 
now possible to contact separately the layers of electrons 
and holes in GaAs, separated by an AlGaAs barrier of 
thickness d as small as 15 nm. 

In this paper, we consider the BCS-BEC crossover in 
electron-hole bilayers when the densities of the electrons 
and holes are varied independently of each other. The ef- 
fect of the density imbalance resembles that of a magnetic 
field in a superconductor (disregarding orbital effects) 
first considered by Sarma [lj]. This analogy was also 
noted in Ref. [lfij , where the influence of an isospin asym- 
metry in nuclear matter was considered for the BCS-BEC 
crossover. The specific system we shall consider as an ex- 



ample is GaAs- AlGaAs. 

Recent experiments with population imbalance in ul- 
tracold trapped Fermi atoms yj| have stimulated a con- 
siderable amount of theoretical work on two-component 
Fermi systems with density imbalance ^3 E3- Em- 
phasis has been placed on the possible occurrence of 
exotic phases in addition to the ordinary BCS pair- 
ing [l3)llij|2fj. However, in these systems the presence of 
a trap and their charge neutrality inhibit the occurrence 
of the exotic phases, and so far only phase separation be- 
tween a superfluid core with equally matched populations 
and an outer normal phase has been detected. 

Electron-hole bilayers may offer a better opportunity 
of observing such exotic phases because the Coulomb re- 
pulsion within each layer acts to suppress phase sepa- 
ration 21]. In addition, we find that the different elec- 
tron and hole effective masses in GaAs, m e and mj,, and 
also the non-local nature of the electron-hole attraction 
both favor the occurrence of exotic phases. These in- 
clude the Sarma phase with one or two Fermi surfaces 
and finite population imbalance, and the Fulde-Ferrel- 
Larkin-Ovchinnikov (FFLO) phase |22f . as well as the 
ordinary BCS pairing with equal populations. The rel- 
ative extension of these phases is quite asymmetric be- 
tween an imbalance with more holes and an imbalance 
with more electrons. We also find that in the superfluid 
phase the separate electron and hole chemical potentials 
display a jump when reversing the population imbalance 
from more electrons to more holes, while in the normal 
phase no jump occurs. Detection of this jump could thus 
serve to identify superfluid character in a system. 

The electron-hole bilayer system is described by the 
grand-canonical Hamiltonian: 
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Here, k, k', and q are two-dimensional wave vectors in 
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the layers, f2 is the quantization volume (surface area), 
ct (ckcr) are the creation (destruction) operators for elec- 
trons (e) and holes (h) distinguished by a = (e, h), the 
^kcr = eko-— /ier are the band dispersions with chemical po- 
tentials \x a for electrons and holes, £kc = k 2 /(2m c ) + E g 
and £h(k) = k 2 /(2mh). The semiconductor band gap 
E g can be reabsorbed in the electron chemical potential. 
Explicit spin quantum numbers are omitted. 

We have carried out our calculations in the zero- 
temperature limit, where a mean-field description of the 
BCS-BEC crossover is appropriate even in two dimen- 
sions 113] ■ We use an unscreened electron- hole attractive 
potential as in Ref. 7], V£ h = -2ne 2 exp (-kd)/(ke), 
where k = |k|, e is the electron charge, and e the 
background dielectric constant. The gap in the super- 
fluid phase makes screening less effective, so introducing 
screening should only cause small quantitative changes. 

In the present calculation we neglect the intra-layer 
Coulomb repulsions V cc and V hh . In the absence of den- 
sity imbalance, V cc and V hh can be readily included in 
a mean- field treatment 0, but when the densities are 
imbalanced, including V cc and V hh in the mean field re- 
quires a detailed knowledge of how overall charge neutral- 
ity is attained. This is in order to avoid divergence of the 
Hartree term. This depends on the specific engineering 
configuration of the device, and therefore in the interests 
of generality we drop V ee and V hh here. When the densi- 
ties are equal we can compare with Ref. , and we have 
verified that omitting V ee and V hh reduces the size of the 
gap by an amount no larger than 30 %. As a separate 
issue, and as we have already noted, intra-layer repulsion 
should stabilize the system against phase separation. For 
this reason, we exclude the possibility of phase separation 
from our discussion of the phase diagram. 

The relevant mean-field equations to be solved for the 
variables /i , /ih, and the (s-wave) gap function Ak are: 

A * = -^E^-k'^t 1 -/^)-/^)], (2) 

k' 
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where f(E) = Q(—E) is the Fermi function at zero tem- 
perature, £ k = + A£ with a = (S + £k)/2. and 
E± = E k ± (5£ k with <5£ k = (££ - fk)/ 2 - In addition, 
i»k = 1 - m 2 . = (1 - £k/_E k )/2. Throughout, we ex- 
press lengths in units of the effective Bohr radius On, 
and energies in units of the effective Rydberg Ry* [24| . 
In two dimensions the average interparticle spacing is 
r s = [TrK + rih)/^- 1 / 2 . 

Figure 1(a) shows the wave- vector dependence of Ak 
for equal densities, layer separation d = 1, and sev- 
eral values of r s . These values of r s and a layer sepa- 
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FIG. 1: (a) Wave-vector dependence of the gap function for 
a = 0, d = 1, and several values of r 3 ; (b) Maximum value 
A max = max{Ak} vs a for d = 1 and several values of r s . 



ration of d = 1 are already experimentally attainable. 
The wave-vector position fc max of the peak in the gap 
function is seen to evolve from a finite value in the BCS 
regime (small r s ) toward zero in the BEC regime (large 
r s ), while the corresponding value A max of the gap func- 
tion attains its maximum for intermediate values of r s . 
This is a generic feature of the density-induced BCS-BEC 
crossover [25^ . 

We have verified for several values of the distance 
1 < d < 4 that the optimal value of r s at which A max 
attains its largest value is located close to the value of 
r s where the average chemical potential /i = (fi c + /ih)/2 
crosses zero. For d = 1 this occurs near r s — 3. This find- 
ing is in line with results for the BCS-BEC crossover with 
a contact potential in three dimensions, namely, that su- 
perfluid properties are more robust in the region located 
between the Fano-Feshbach resonance and the vanish- 
ing of the chemical potential [2(|. The narrow region 
lying between the optimal value of r s and the value cor- 
responding to /i = thus identifies the middle of the 
crossover region between the BCS and BEC regimes for 
the electron-hole bilayer. 

The effect of the density imbalance a = (n e —ri} l )/(ri e + 
7ih) on A max is shown in Fig. 1(b). We see that the den- 
sity imbalance acts to reduce the magnitude of the energy 
gap and that it has different effects on the two sides of 
the crossover. In the BCS regime, the mismatch of the 
Fermi surfaces for electrons and holes strongly affects the 



3 



1 

0.5 
a 
-0.5 
-1 



N 


/ 1 

' 




52/ I 




, S1 


___---fFI_C 




N \ 


si \ 







0.5 1 



2 4 



8 16 



FIG. 2: Zero-temperature phase diagram for d = 1, showing 
the stability domains of the different phases: Normal (N), 
Sarma (Sf and S2), and FFLO. The dashed line corresponds 
to the curve p — 0, and the dash-dotted line separates the SI 
and S2 phase. Note the logarithmic scale for r s . 



superfluid properties of the system. The superfluid prop- 
erties are lost when the Fermi energies mismatch becomes 
large compared with the value of A max for equal densi- 
ties. For r s small this occurs for a small value of a. In 
the BEC regime, the superfluid properties are less sen- 
sitive to density imbalance, whose main effect is then to 
reduce the number of electron-hole pairs. 

Figure 2 shows the zero-temperature phase diagram for 
d = 1. We can identify various phases using Ak, deter- 
mined from Eq. (2), and the superfluid (mass) density p s . 
Within mean-field theory and in the zero-temperature 
limit, p s is given by: 
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Here, kj is the j-th zero of E% — with A = ( 



-). [For 

positive (negative) density imbalance only E£ (E^ ) has 
zeros, while no zero occurs for a = 0.] 

The normal phase (N) corresponds to the trivial solu- 
tion Ak = 0. The Sarma phases corresponds to nonva- 
nishing Ak when a ^ and positive superfluid density 
p s . The SI and S2 denote the Sarma phases for one 
and two Fermi surfaces, respectively. There will be one 
zero of Eh (j = 1) for the SI phase (one Fermi surface), 
and two zeros (j = 1,2) for the S2 phase (two Fermi 
surfaces). Sometimes, the Sarma S2 phase is called the 
"breached-pair" phase after Ref. 19]. A negative value of 
p s in Eq. (5) indicates that the Sarma phase is unstable 
toward a phase with a spontaneous ly g enerated super- 
fluid current, which we associate |27ll28|| with the FFLO 
phase. We have verified that the Sarma phase, whenever 
it exists, is always lower in energy with respect to the 
normal phase. (We recall in this respect that our cal- 
culation is at fixed density imbalance, while the original 
Sarma calculation was at fixed chemical potentials.) 



The most prominent feature of Fig. 2 is the marked 
dependence of the phase diagram on the sign of a [29j |. 
In particular, while the boundary of the normal phase 
does not depend appreciably on the sign of a, the region 
of stability of the Sarma phase with respect to the FFLO 
phase depends dramatically on the sign. For a < 0, the 
phase diagram is dominated by the FFLO phase, with 
the Sarma (SI) phase being confined to the extreme BEC 
region, while for a > 0, the FFLO phase is compressed 
into the region of small r s ■ 

Such an asymmetry can be understood in terms of the 
relevant dispersion E£ (E^) for positive (negative) a. 
Due to the mass difference, the term 2<5£k which is added 
(subtracted) to Ek makes the dispersion E£ (££ ) steeper 
(flatter) with respect to the case of equal masses. As it 
is clear from the last term in Eq. (J5J, a flatter disper- 
sion will make the superfluid density more negative, thus 
replacing the Sarma phase with the FFLO phase. The 
opposite occurs for positive a, when the relevant disper- 
sion becomes steeper. 

In this part of the phase diagram there is some room 
even for the S2 Sarma phase. So far, this phase was found 
in the literature to be very fragile, being invariably re- 
placed by phase separation or by the development of a 
FFLO phase [ni|22l|3fj. Here, the concurrence of several 
favorable factors, the mass difference, the wave- vector 
dependence of the gap, and the intralayer Coulomb re- 
pulsion, all serve to stabilize the S2 Sarma phase in an 
appreciable region of the phase diagram. 

We note that the largest number of phases occurs for 
intermediate values of r s (say, r s — 1.5 — 3), correspond- 
ing to the smooth crossover region between the BCS and 
BEC limit for a = 0. For non-zero a, several transition 
lines can be crossed in this region by varying a or r s . 
This is thus the most fertile region to be explored experi- 
mentally, also because we recall in the same region A max 
is largest and the superfluidity most robust. 

In this respect, the behavior of the separate chemical 
potentials p c and ph vs a can serve to reveal the appear- 
ance of a finite value for the gap A. This is because p e 
and /Lth must show a jump across a = in order to sus- 
tain a finite density imbalance in the superfluid phase. 
This behavior is shown in Fig. 3 for d = 1 and r s = 3, 
when p e (a = + ) — p c (a = 0~) ~ 2 A max . The possible 
instability toward the FFLO phase we discussed above 
should not affect the occurrence of this jump. Physical 
quantities obtained for the Sarma and the FFLO phases 
should, in fact, merge continuously when a approaches 
zero, owing to the corresponding vanishing of the modu- 
lation wave vector associated with the FFLO phase [13] . 

It is useful to discuss the sensitivity of our results to 
the parameters used in the calculations. Increasing the 
distance d will shift the location of the intermediate re- 
gion where most phases are seen, to larger values of r s 
(for d = 2, e.g., the shift is about a 30%), but it does not 
alter the shape of the phase diagram. Reducing the mass 
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FIG. 3: Zero-temperature electron (fi e ) and hole (fj,^) chemi- 
cal potentials vs a for d = 1 and r s = 3. The average chemical 
potential fi is also shown. 



difference has the effect of contracting the FFLO phase 
for negative a and expanding it for positive a. This is as 
expected from the arguments discussed above. However, 
even for equal masses the phase diagram changes only 
quantitatively, allowing even in this case some space for 
the S2 Sarma phase. Finite temperature assists in stabi- 
lizing the Sarma phases with respect to the FFLO: the 
superfluid density is, in fact, quite sensitive to tempera- 
ture in the presence of density imbalance. We have found, 
however, that temperatures below 1 K do not appreciably 
alter the phase diagram in the most interesting interme- 
diate r s region. Finally, a reduction of the FFLO phase 
due to the intra-layer Coulomb repulsion may occur but 
this reduction should be limited to the region of large 
r s , where the spatial modulation of the gap parameter is 
expected to be accompanied by a density modulation. 

In conclusion, we have shown that electron-hole bilay- 
ers are promising candidates for revealing a variety of su- 
perfluid phases while traversing the BCS-BEC crossover. 
Exotic phases such as the FFLO and the S2 Sarma phases 
which have been so far elusive to experimental detection, 
should be notably robust in electron-hole bilayer systems. 
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